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ABSTRACT 

We report the identification of the double degenerate system NLTT 16249 that comprises a normal, 
hydrogen-rich (DA) white dwarf and a peculiar, carbon-polluted white dwarf (DQ) showing photo- 
spheric traces of nitrogen. We disentangled the observed spectra and constrained the properties of 
both stellar components. In the evolutionary scenario commonly applied to the sequence of DQ white 
dwarfs, both carbon and nitrogen would be dredged up from the core. The C/N abundance ratio 
(w 50) in the atmosphere of this unique DQ white dwarf suggests the presence of unprocessed mate- 
rial ( 14 N) in the core or in the envelope. Helium burning in the DQ progenitor may have terminated 
early on the red-giant branch after a mass-ejection event leaving unprocessed material in the core 
although current mass estimates do not favor the presence of a low-mass helium core. Alternatively, 
some nitrogen in the envelope may have survived an abridged helium-core burning phase prior to 
climbing the asymptotic giant-branch. Based on available data, we estimate a relatively short or- 
bital period (P < 13 hrs) and on-going spectroscopic observations will help determine precise orbital 
parameters. 

Subject headings: binaries: close — stars: individual (NLTT 16249) — white dwarfs 



1. INTRODUCTION 

The presence of a large concentration of carbon, de- 
tected via C2 molecular bands, in the otherwise helium- 
rich atmosphere of many old (> 10 9 years) white 
dwarfs has been viewed for many years (e.g., Vauclair 
& Fontaine 1979; Koester et al. 1982) as evidence of the 
carbon-rich nature of their cores: The material diffuses 
upward, away from the core, and is dredged-up to the 
surface by the deep helium convection zone developing 
in cool white dwarfs (Fontaine et al. 1984; Pelletier et 
al. 1986; MacDonald et al. 1998). These objects, collec- 
tively known as DQ white dwarfs, form a homogeneous 
sequence of stars with surface temperatures ranging from 
~ 11, 000 to <~ 5, 000 K and with temperature-correlated 
surface carbon abundances between 10 -2 and 10 -7 rela- 
tive to helium by number (Dufour et al. 2005; Koester & 
Knist 2006). 

Evolutionary models (e.g., Schaller et al. 1992) show 
that stars of low to intermediate masses, i.e., below 8 to 
10 Mq, end their nuclear burning history with a core 
composed primarily of carbon and oxygen or heavier 
elements (Ne, Mg) in the upper mass range (Garcia- 
Berro et al. 1997). In the past, spectroscopic searches 
for oxygen (ultraviolet CO bands) or nitrogen (violet 
CN bands) have yielded abundance limits 0/C« 0.1 — 1 
increasing with temperature (Koester et al. 1982) or 
0/C,N/Cw 10- 2 -10~ 3 (Wegner & Yackovich 1984). Al- 
though oxygen is a primary product of the triple-alpha 
(3-a) nuclear burning process, the actual amount of oxy- 
gen dredged-up along with carbon from the interior de- 
pends on its abundance profile in the core (MacDonald et 

'Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile under 
programme ID 086.D-0562. 

2 Visiting Astronomer, Kitt Peak National Observatory, Na- 
tional Optical Astronomy Observatory, which is operated by the 
Association of Universities for Research in Astronomy (AURA) un- 
der cooperative agreement with the National Science Foundation. 



al. 1998). On the other hand, nitrogen should be fully ex- 
hausted along with helium by the 3-a process and is not 
expected in the atmosphere of white dwarfs with a fully- 
evolved core. Recently, Gansicke et al. (2010) reported 
the discovery of two oxygen-rich white dwarfs alongside 
the carbon-rich DQ sequence. The measured O/C abun- 
dance ratios (log O/C « 0.6 — 1.8) imply core composi- 
tions dominated by oxygen (80 to 99% by mass). The 
stars were suspected by Gansicke et al. (2010) to be the 
descendents of the most massive stars to avoid core col- 
lapse (9-lOM ; Garcia-Berro et al. 1997). 

Clearly, the surface composition reflects the white 
dwarf core composition which, in turn, is determined by 
the evolutionary path followed by the progenitor. In this 
context, we present in Section 2 spectroscopic observa- 
tions of the peculiar object NLTT 16249 that show a 
strong Balmer line series along with C2 Swan and CN 
violet band absorptions. The observed radial velocity 
variations and a successful spectral decomposition and 
model atmosphere analysis (Sections 3.1-3.3) show that 
the object is in fact a close double degenerate system 
comprising a normal DA white dwarf and a peculiar 
DQ of similar optical luminosities. We examine possi- 
ble evolutionary scenarios in Section 3.4, and we propose 
that the photospheric carbon to nitrogen abundance ra- 
tio (C/Nrj 50) indicates that the core material was left 
partially unprocessed or that burning ceased before de- 
stroying all nitrogen in the envelope. Finally, we sum- 
marize in Section 4. 

2. OBSERVATIONS AND CLASSIFICATION 

We first observed NLTT 16249 on UT 2010 Mar 26 
with the R-C spectrograph attached to the 4m telescope 
at Kitt Peak National Observatory. We obtained three 
consecutive exposures of 900 s each using the KPC-10A 
grating with the WG360 order blocking filter. The spec- 
tra cover a wavelength range from 3660 to 6790 A with 
a dispersion of 2.8 A per pixel. This low dispersion spec- 
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Fig. 1 — Spectrum of NLTT 16249 obtained with the UVB arm 
of X-shooter. The main spectral features (Balmer lines and C2 or 
CN molecular bands) are marked with dashed lines and labelled 
accordingly (see text). 



trum (A/AA w 1000) revealed a spectroscopic combina- 
tion that would be unusual in any star. The hydrogen 
Balmer line series are blended with C 2 Swan bands (the 
d 3 TL g — a 3 H u system; see Tanabashi et al. 2007) and what 
appears to be the Aw = sequence of the CN violet band 
near 3880 A (the B 2 Y,+ - X 2 Y>+ system; see Ram et al. 
2006). 

Next, we obtained a single echelle spectrum on UT 
2010 Nov 7 (2400 s exposure) with the X-shooter spectro- 
graph (Vernet et al. 2011) attached to the UT2 (Kueyen) 
at Paranal Observatory. The spectrum covers a wave- 
length range from 3000 A to 2.5 ^m on three sepa- 
rate arms with a resolving power A/AA « 9100 for the 
UVB arm, 8800 for the VIS, and 6200 for the NIR. This 
intermediate dispersion spectrum confirmed our earlier 
findings and also revealed the presence of additional se- 
quences from the B — X CN system. 

The star could be identified as a DQAZ, i.e., a helium- 
rich white dwarf with traces of carbon and hydrogen 
along with additional elements such as nitrogen, al- 
though the strength of the hydrogen lines is incompat- 
ible with this interpretation. More likely, NLTT 16249 
is a DA+DQ(N) double degenerate system comprised of 
a normal hydrogen-rich (NLTT 16249B) and a nitrogen- 
contaminated DQ white dwarf (NLTT 16249A). We pro- 
pose a non-standard classification for the DQ(N) white 
dwarf owing to the unforeseen presence of nitrogen (for 
the current classification scheme sec Sion et al. 1983). 
The diagnostics employed to confirm the second possibil- 
ity are, first, a radial velocity study and, next, a spectral 
decomposition into two components. 

Figure 1 shows the X-shooter spectrum of NLTT 16249 
(UVB arm): The main molecular bands are marked at 
selected wavelengths (Wallace 1962) following standard 
spectroscopic notations along with the Balmer lines. 

The relative velocity change between the two compo- 
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nents was readily measurable by using spectra obtained 
at the two available epoches, although individual ve- 
locity changes may suffer from large uncertainties on 
the zero-point of the low-dispersion KPNO wavelength 
scale. The relative velocity change between the two 
components is free of systematic errors and allowed us 
to estimate i^DQ + ^da ^ 320 ± 30 kms -1 . On the 
other hand, the velocity shift of the DQ component was 
u i.dq~ w 2,dq+C = 120±30kms _1 , where the subscripts 
"1" and "2" refer to the KPNO and VLT spectra, respec- 
tively, and C is the zero-point error. Similarly, for the 
DA component t;i ; DA—V2,DA+C = — 200±30kms~ 1 , i.e., 
in opposite direction to the DQ component. The ampli- 
tude ratio would imply that the DA white dwarf is more 
massive than its companion with Mdq/Mda = l-7^ 6 , 
but because of potential systematic errors we will simply 
assume that Mdq ~ -Mda- 

The binary has a significant proper motion of 
/i Q cos S = 2±6, us = -140±6 mas yr" 1 (Salim & Gould 
2003) in a relatively crowded field. We collected Johnson 
photometry (V = 15.77, B-V = 0.24, U - B = -0.66 
with an assumed error of 0.02 on the colors) from Eggen 
(1968) and 2MASS photometry (J = 14.87 ± 0.04) from 
Skrutskie et al. (2006). A close examination of the VLT 
acquisition image (epoch 2010) shows a fainter crowded 
star to the NW. Taking into account the proper motion 
of NLTT 16249, the stars would have nearly overlapped 
circa 1990. The 2MASS-J frame (epoch 1997) clearly 
shows a SE-NW elongation. Also, the DSS1 (epoch 
1955) red and blue plates show that the crowded star 
is markedly red so that it probably contaminated the 
2MASS-J photometric measurement. With the crowded 
star at a larger separation, the UBV photometry ob- 
tained by Eggen (1968) circa 1965 was probably not 
significantly contaminated. In order to extend avail- 
able broadband indices further in the red, we folded 
the X-shooter UVB and VIS spectra with VRI band- 
passes (Bessell 1990) and measured V — R = 0.15 and 
V - 7 = 0.37. 

3. ANALYSIS 
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Fig. 2. — Johnson-Cousin color indices B — V versus V — I (left panel) and U — B versus B — V (right panel) for pure-hydrogen 
(DA) model grids (full lines) and DQ model grids (dashed lines). The DA model grids are shown with varying surface gravity (labelled 
at logg = 7.5,8.0, and 8.5) and temperature (shown from right to left at 7, 8, 9, 10, 12, 14, 16xl0 3 K). The DQ model grids are shown 
with varying carbon and nitrogen abundance (labelled with the pair log C/He, log N/He) and temperature (shown from right to left at 7, 
7.5, 8xl0 3 K) at log<? = 8. On both panels the observed colors arc shown with error bars, and the best-fit model colors are shown, from 
upper-left to lower-right, with grey squares and lines at (log C/He, log N/He) = (-4,-5.7), (-4.5,-6.2), (-5.0,-6.7), (-5.5,-7.2). 



The components were analyzed conjointly with the to- 
tal flux measured at earth given by 

U = ^ (-^DQ Hv,T>Q + RdA #^,Da) (1) 

where D is the distance to the binary, R their respective 
radii, and H v their respective model Eddington fluxes. 
The individual radii, hence masses, are constrained by 
the surface gravity of each model grid point using mass- 
radius relations from Benvenuto & Althaus (1999). The 
best-fit parameters (T cffjDQ , logg DQ , T eff;D A, log g da, 
and DQ abundances C/He, N/He) are obtained using \ 2 
minimization techniques. We also computed Johnson- 
Cousin UBVRI magnitudes and color indices of the com- 
posite spectra to constrain the best-fit parameters. Be- 
cause of the absence of CO spectroscopic signatures we 
set the oxygen abundance to zero in the DQ models. 
Moreover, the absence of Ca H&K and CH G-band also 
suggests very low metallicity in cither star and the ab- 
sence of hydrogen in the DQ white dwarf. With the re- 
vealing exception of nitrogen the DQ component appears 
normal. 

3.1. Model atmospheres 

The DA model atmospheres are described in Kawka & 
Vennes (2006, 2011). The new DQ model atmospheres 
are in dual convective/radiative equilibrium as well as 
in local thermodynamic equilibrium. We included all 
relevant species (He, He+, C, C+, N, N+, O, 0+) and 
molecules (C 2 , CN, CO) in the charge conservation and 
abundance equations: This non-linear system of equa- 
tions was solved using the secant method before each 
model iteration. We employed the molecular partition 
functions of Sauval & Tatum (1984). We also estimated 
the contributions of C^~ and concluded that virtually all 
electrons are contributed by the ionization of helium and 
carbon. 



The models add all relevant opacities, including ul- 
traviolet C I lines, and the C2 Swan bands (Av = 
-1, 0, +1, +2) and the violet CN bands (Av = -1, 0, +1) 
using the "just-overlapping line approximation" (Golden 
1967; Zeidler-K.T. & Koester 1982). We also added the 
Hc~ free-free and the C~ bound-free and free-free opac- 
ities as well as helium Rayleigh scattering. We updated 
the Zeidler-K.T. (1987) CN molecular opacities with os- 
cillator strengths from Bauschlicher et al. (1988) and 
molecular constants tabulated in Reddy et al. (2003) in 
good agreement with Ram et al. (2006). Table 1 lists 
adopted bandhead wavelengths and gf values for the first 
three vibrational transitions of the C2 and CN systems. 

Figure 2 shows relevant photometric properties of the 
model grids. In particular, optical DQ colors are sensitive 
to carbon (and nitrogen) abundance because of varying 
depths of Swan and CN violet bands. This effect influ- 
ences the composite colors and the selection of best-fit 
composite models. We found that the DQ spectral en- 
ergy distribution is only mildly affected by surface grav- 
ity variations (see Dufour et al. 2005). The DQ synthetic 
colors are markedly bluer than those of a DA model at 
the same temperature. 

3.2. Properties of the components 

First, we proceed with the spectral decomposi- 
tion. We mapped the minimum x 2 values in the 
(log <7dq, log C/He) plane corresponding to the best-fit 
DA parameters (T^da, log <7da) and DQ temperature 
(Tefl^DQ)- We varied log C/He from —5.5 to —4.0 and 
logSDQ from 7.5 to 8.75. In all spectral decomposi- 
tions the nitrogen abundance scaled with the carbon 
abundance with a ratio logC/N m 1.7 that we subse- 
quently held fixed at this value. The calculations show 
that a family of solutions exists along an axis in the 
(7eff,DQ, log C/He) plane where higher DQ temperatures 
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Fig. 3. — (Left panel) Model decomposition of the observed spectrum (degraded to a resolution of 2.4A for clarity). The model flux 
contributions of the DA and DQ components are shown with dashed and dotted lines, respectively. The total DA+DQ model flux (full 
line) is compared to the observed flux (grey line). (Right panel) Composite model spectrum fit to the X-shooter spectrum (full resolution) 
near the Balmcr scries limit and the CN violet band (Av = 0). 



are compensated with higher carbon abundances. This 
degeneracy of the solutions may be lifted by comparing 
the predicted colors of the model composites with ob- 
served colors. On the other hand, the minimum \ 2 val- 
ues in the (logc/DQ,logC/He) plane are steadily found 
for log g da ~ 8.2 — 8.5. Interestingly, the best-fit surface 
gravities (hence radii) always indicated similar binary 
component luminosities, so that the DQ gravity corre- 
lates with the DA gravity and was found in the range 
log5D Q «8.0-8.5. 

The predicted color indices were compared to observed 
indices to constrain possible solutions. Figure 2 shows 
composite colors corresponding to the best-fit models 
with logC/He varying from —5.5 to —4.0. As noted 
above, higher carbon abundances are compensated with 
high DQ temperatures. Both sets of indices favor a lower 
carbon abundance 

logC/He = -5.0 ±0.3, logN/He = -6.7 ±0.3, 

and we conclude that the DA model parameters are 

Toff, da - 8330 ± 200 K, \ogg DA = 8.37 ± 0.15, 

and, correlating with the DA parameters, the DQ pa- 
rameters are 

T cS ,bq = 7770 ± 200 if, log ffDQ = 8.27 ± 0.22. 

The total system mass varies between 1.33, in the lower 
gravity range, to 1.83M , in the higher gravity range. 
The distance calculated using the V — My modulus varies 
from ~34 pc for the high system mass to ~46 pc for 
the low system mass. A parallax measurement would 
help confirm and narrow down the range of spectroscopic 
solutions. 

Assuming a system mass in the range Mda + ^dq = 
1.3 — 1.8Mq and a minimum velocity amplitude if da + 
Kdq > 320 kms -1 , the orbital period is constrained to 

M D a + M D q .3 

sin i < 9.2 — 12.8 hrs. 



P = 2irG 



(K UA + if DQ ) E 



Owing to the large velocity amplitude and relatively 
short period, the orbital parameters (P, e, if da, ^dq) 
should be easy to determine. 

Figure 3 shows the best-fit composite model spectrum 
to the X-shooter spectrum using parameters listed above. 
The Swan Aw = +1 bands appear too strong relative to 
Aw = bands, particularly in the extended wing. For- 
tunately, the abundance offset between these two bands 
does not exceed s=s 0.2 dex. 

3.3. The dominant carbon isotope 

The isotopic shift of the Av = +1 band of 13 C 12 C rela- 
tive to that of C 2 is 7.5 A (Wegner 1984). By attributing 
the existing band to C 2 we conclude that the absence 
of red-shifted absorptions implies C2/ 13 C 12 C > 10, i.e., 
12 C/ 13 C>10. 

3.4. Evolutionary scenarios 

Hydrogen-deficient white dwarfs (DB and DQ white 
dwarfs, see Sion et al. 1983) are possibly the product of 
a late thermal pulse. This process leads to the formation 
of a hydrogen-deficient surface that is also enriched in 
carbon and oxygen (Herwig et al. 1999, and references 
therein) . After reaching the white dwarf cooling sequence 
these objects develop a nearly pure-helium surface and 
are recognized as DB white dwarfs. While cooling further 
they develop a deep helium-dominated convection zone 
dredging-up carbon from the core (Fontaine et al. 1984; 
Pelletier et al. 1986; MacDonald et al. 1998). 

The material mixed in the helium-dominated convec- 
tion zone of the DQ white dwarf NLTT 16249A is en- 
riched with carbon and nitrogen. The diffusion time- 
scales at the bottom of the helium convection zone of 
a <~ 8000 K white dwarf are relatively short (w 10 3 
yr; Koester 2009). Therefore, a steady-state abun- 
dance of nitrogen and carbon is achieved only if the 
carbon/nitrogen-rich material is supplied to the envelope 
from the core. This scenario, commonly applied to ordi- 
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nary carbon-rich DQ white dwarfs, should also apply to 
the peculiar DQ NLTT 16249A, but the source of nitro- 
gen remains elusive. 

A copious amount of nitrogen is produced in 
intermediate-mass stars between pulses on the AGB 
when the convective envelope reaches the hydrogen- 
burning shell (Lattanzio et al. 1996). This production 
is deemed responsible for nitrogen enrichment in plan- 
etary nebulae (Karakas et al. 2009), but what fraction 
of it would remain within the helium envelope and con- 
tribute to the white dwarf core composition is an open 
question. 

Nitrogen leftover material is present in the early-stages 
of helium burning on the horizontal-branch (HB), but it 
is quickly destroyed along with helium leading to the 
formation of a carbon/oxygen-rich core. For example, a 
carbon-to-nitrogen ratio of « 50 is briefly achieved in the 
core of 2 Mq HB star after ~ 110 Myr, but it rapidly de- 
clines afterward (Schaller et al. 1992). At this particular 
stage the core composition is 64.4% He, 27.6% C, 0.6% N, 
and 7.4% O by mass, but by the end of the helium burn- 
ing phase the remainder consists mostly of oxygen (83.6% 
by mass) and carbon (12.7%) without a trace of nitrogen. 
Prior to helium ignition and on the RGB the carbon- 
to-nitrogen number ratio in the core is only w 1/100. 
At higher masses and up to 9 Mq nitrogen is destroyed 
even more rapidly during the core helium-burning phase, 
and on the RGB the carbon-to-nitrogen number ratio in 
the core never rises above ~ 1/40. Should burning be 
terminated while on the RGB, nitrogen would dominate 
carbon in a helium-rich core. 

We found that the DQ NLTT 16249A is part of a 
close double degenerate system. Binary evolutionary sce- 
narios (e.g., Nelemans et al. 2001) show that, normally, 
the mass of helium-core white dwarfs correlate with the 
double-degenerate orbital period and almost never ex- 
ceed 0.4 Mq. Our spectral decomposition suggests a DQ 
mass larger than w 0.6 Mq. Applying the Schonberg- 
Chandrasckhar mass limit to the case of a helium core 
and helium/hydrogen envelope, the core mass on the 
RBG could reach M c /M = 0.08 prior to collapse and 
ignition, i.e., M c w 0.6 Mq in the case of a M = 8 Mq 
star. In order to conceal a large helium-rich core the 
DQ white dwarf in NLTT 16249 would have to be the 
product of such exceptional circumstances. 

Based on our spectral decomposition we find that the 
stars have comparable luminosities with similar masses 
and cooling ages. In the lower mass range the cooling 
ages are t dq = 1.3 Gyr and £da = 1-5 Gyr so that 
the DA may have preceded the DQ by < 200 Myr, but 
in the higher mass range the situation is reversed with 
t dq = 2.6 Gyr and £da = 2.4 Gyr. In either case, 
the progenitors would follow each other off the main- 
sequence within 200 Myr or less. The maximum age dif- 
ferential could apply to initial masses of 3 and 4 Mq , or 
in a lower mass range to initial masses of 2.0 and 2.2 Mq 



(see Schaller et al. 1992). Narrowing the age differential 
would accommodate even larger initial masses and we 
cannot exclude the possibility of a M — 8 Mq progeni- 
tor for the DQ white dwarf although it is less likely than 
a 2-4 Mq progenitor. Moreover, a high-mass progenitor 
for the 0.6-0.9 Mq DA white dwarf is also unlikely (Wei- 
demann 2000). Therefore, NLTT 16249 is most probably 
a CO+CO system with the DA forming first, and the DQ 
forming next during a common-envelope phase; the pro- 
genitor masses would be in the 2-4 Mq range. 

The question of the source of nitrogen remains. So 
far, we considered the composition at the center of the 
star as representative of the core. As pointed out by 
MacDonald et al. (1998) in the case of oxygen, the ma- 
terial actually dredged-up is not located at the core cen- 
ter but, instead, below the former helium-burning shell. 
Therefore, we propose that helium burning terminated 
prematurely during a common-envelope phase leaving an 
abundance ratio C/Nss 50 off-center and below the for- 
mer helium-burning shell in NLTT 16249 A. In conclu- 
sion, the presence of nitrogen in the atmosphere of DQ 
white dwarfs would help distinguish normal DQ white 
dwarfs that evolved in isolation and DQ white dwarfs 
that are the product of binary evolution. 

4. SUMMARY 

We report the discovery of a peculiar DQ white dwarf 
in a close degenerate system (P < 13 hrs). The pres- 
ence of nitrogen in the atmosphere of the DQ also sug- 
gests the presence of unprocessed material in the core of 
the star, most probably off-center and below the extinct 
helium-burning shell. Helium burning may have been in- 
terrupted following a common-envelope event leaving a 
nitrogen-rich layer above a normal carbon/oxygen core. 
Our spectral decomposition suggests a system mass in 
the 1.3 to 1.8 Mq range with cooling ages between 1.3 
and 2.6 Gyr depending on the mass. The binary will 
merge within 9 Gyr for a system mass of 1.3 Mq, or 13 
Gyr for 1.8 M© following Ritter (1986). New radial ve- 
locity measurements will be used to set a precise binary 
mass ratio and confirm the evolutionary prospects of the 
binary NLTT 16249. 
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